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A  recent  theoretical  study  of  double  Langmuir  probes  led  to  development  of  improved 
analytical  techniques  that  account  for  probe  electrode  sheath  expansion  in  a  low- 
temperature  plasma  environment.  The  proposed  technique  is  based  on  analytical  curve 
fitting  to  Laframboise’s  results,  which  enables  double  probe  measurement  of  electron 
temperature  and  plasma  density  over  a  wide  range  of  densities  without  advance  knowledge 
of  the  probe  radius  relative  to  Debye  length.  In  this  investigation,  the  method  was  evaluated 
by  comparing  plasma  properties  from  three  double  Langmuir  probes  of  varying  electrode 
diameter.  Far-field  measurements  of  a  Hall  thruster  plume  were  conducted  with  each  probe 
at  two  facility  background  pressures.  Plasma  characteristics  were  calculated  using  the 
proposed  analytical  technique,  and  compared  to  orbital  motion  limited  and  thin-sheath 
analyses  in  regions  of  the  plume  where  these  methods  are  valid.  Results  revealed  equivalent 
plasma  properties  from  each  double  probe  geometry  when  using  the  proposed  double  probe 
analysis  methodology,  indicating  the  electrode  sheath  expansion  is  accounted  for.  Calculated 
plasma  densities  and  electron  temperatures  were  consistent  with  standard  analysis 
techniques.  These  findings  support  the  proposed  double  probe  analysis  method  for 
examination  of  plasma  properties  over  a  wide  span  of  the  far-field  plume  using  a  fixed 
Langmuir  probe  design  and  a  single  analysis  technique. 

I.  Introduction 

THE  growing  use  of  electric  propulsion  for  efficient,  in-space  maneuverability  drives  evaluation  of  electric 
propulsion  plumes  in  ground-based  facilities  [1,2,3].  These  investigations  require  measurement  over  increasing 
range  of  plasma  properties  for  laboratory  propulsion  systems  spanning  orders  magnitude  in  power,  from  Watts  to 
over  lOOkW  [1,4].  In  addition,  properties  of  the  internal  thruster  plasma  and  near- field  regions  are  orders  of 
magnitude  higher  than  the  extended  far-field  plume.  Langmuir  probes,  an  electrostatic  diagnostic  developed  by 
Irving  Langmuir  in  1924  [5],  are  widely  used  to  examine  plasma  properties  in  laboratory  investigations  and  on-orbit 
environments.  Although  the  electrostatic  theory  of  Langmuir  probe  operation  is  straightforward,  the  implementation 
and  analysis  is  complicated  by  several  factors,  including  probe  sheath  effects  that  influence  electrode(s)  current 
collection  and  perturbations  to  the  local  plasma.  These  complications  necessitate  the  probes  are  designed  for  a 
specific,  narrow  range  of  plasma  properties  and  require  analysis  techniques  that  are  based  on  assumptions  about  the 
probe  sheath  interactions  with  the  plasma  environment  under  examination  [6]. 
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Langmuir  probes  are  typically  categorized  according  to  geometry  and  number  of  electrodes.  Probe  geometries 
range  from  spherical,  planar,  and  cylindrical  in  nature  with  one  or  multiple  electrodes.  Multiple  electrodes  typically 
decrease  plasma  perturbation  due  to  lower  bias  voltages  at  the  expense  of  increased  probe  size  and  decreased  spatial 
resolution.  The  double  Langmuir  probe  configuration,  with  two  electrodes,  floats  with  the  plasma  and  electrodes  are 
swept  relative  to  the  floating  potential.  This  has  several  advantages  over  the  single  probe  geometry;  it  eliminates  the 
need  for  a  fixed  reference  potential,  limits  biased  electrode  collected  current  to  the  ion  saturation  current,  minimizes 
local  plasma  perturbations  with  a  net  zero  current  to  the  electrode  pair,  and  reduces  the  electrode  voltage  sweep. 
While  the  triple  probe  does  not  require  a  voltage  sweep,  the  separation  between  electrodes  reduces  spatial  resolution. 

Recent  theoretical  studies  of  double  Langmuir  probe  characteristics  in  low-temperature  plasma  led  to 
development  of  a  new  analysis  technique  that  employs  analytical  fits  to  Laframboise’s  numerical  results  and  circuit 
analysis  to  formulate  relationships  between  applied  electrode  voltage,  current  collection,  and  the  local  plasma 
properties  [7].  These  coupled  equations  account  for  probe  sheath  effects,  such  that  it  is  applicable  over  a  wide  range 
of  electron  Debye  length  relative  to  probe  dimensions  and  spans  the  expands  of  measureable  plasma  properties 
compared  to  conventional  analyses. 

In  this  paper,  the  new  double  Langmuir  probe  analysis  technique  [7]  is  assessed  with  measurements  of  a  low- 
power  Hall  thruster  plume  using  three  cylindrical  double  Langmuir  probes  of  varying  scales.  This  systematic  plume 
mapping  encompassed  a  wide  range  of  plasma  properties.  Plasma  parameters  were  calculated  using  the  proposed 
technique  and  compared  to  conventional  methods,  such  as  orbital  motion  limited  (OML)  and  thin  sheath  analysis. 
Plume  mapping  was  conducted  from  20  to  150  thruster  diameters  downstream  at  two  facility  background  pressures. 


II.  Probe  Theory 


Analyses  of  Langmuir  probe  current- voltage  (I-V)  characteristics  are  dependent  on  the  length  scales  of  electron 
Debye  length  and  probe  geometry.  The  size  of  the  plasma  sheath  surrounding  the  probe  is  proportional  to  the 
electron  Debye  length  XD,  as  defined  in  Eq.  (1),  where  Te  is  electron  temperature,  ne  is  electron  number  density,  kB  is 
the  Boltzmann  constant,  and  e  is  the  electron  charge.  The  probe  is  in  a  collisionless,  “thick”  sheath  regime  when  the 
ratio  of  probe  radius  rp  to  Debye  length  is  rf/AD<3,  often  referred  to  as  orbital  motion  limited  (OML).  A  “thin” 
sheath  analysis  is  applicable  for  rp/XD  >10. 
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The  transitional  regime  between  the  “thick”  and  “thin”  sheath  analyses  where  3<r//2Z)<10  may  be  analyzed 
using  an  iterative  approach  with  Laframboise’s  comprehensive  compilation  of  collected  ion  current  with  respect  to 
electron  Debye  length  for  cylindrical  or  spherical  probe  geometry  [8].  Several  analytical  fits  to  Laframboise’s 
results  have  been  proposed,  including  studies  by  Steinbruchel,  Narasimhan,  and  Karamcheti  [9,10,11]  that  led  to 
development  of  analytical  parameterization  techniques  of  Laframboise’s  results  that  were  successfully  demonstrated 
with  single  Langmuir  probe  geometries  for  cold  ions  (T/Te«  1)  over  a  wide  span  of  rp/XD,  where  Tt  is  the  ion 
temperature.  The  parameterization  technique  has  been  reformulated  for  analysis  of  double  Langmuir  probes  [7],  and 
is  valid  for  all  rP/XD  if  the  electron  distribution  function  is  Maxwellian.  Experimental  assessment  of  the 
parameterization  approach  for  double  probe  analysis  is  the  aim  of  this  investigation.  An  overview  of  the  technique 
is  described  below,  and  additional  details  are  provided  elsewhere  [7]. 

A  floating  double  Langmuir  probe  circuit  with  ion  and  electron  currents  is  represented  in  Fig.  1.  The  positive 
sense  of  ion  current  and  electron  current  is  described  by  the  notation  IJ  +  and  Ij>e,  where  j  indicates  electrode  1  or  2. 
Kirchoff  s  laws  are  expressed  in  Eqs.  (2)  and  (3),  where  the  applied  electrode  voltage  Vj  is  relative  to  the  local 
plasma  floating  potential  and  the  magnitude  of  the  probe  current  Ip  cannot  exceed  the  ion  saturation  current.  Both 
electrodes  are  electron  repelling  ( Vj ,  V2<0)  since  the  ion  saturation  current  limit  is  always  much  smaller  than  the 
electron  saturation  current.  One  electrode  is  biased  slightly  below  plasma  floating  potential  and  the  other  slightly 
above.  These  characteristics  enable  formulation  of  electron  currents  Ij>e  in  Eq.  (4),  where  Aj  is  the  surface  area  of  a 
single  electrode,  e  is  the  magnitude  of  the  electron  charge,  n0  is  the  local  number  density  of  the  undisturbed  plasma, 
me  is  the  electron  mass,  f  is  the  local  electron  temperature  in  electron  volts  (i.e.  kBTe/e ),  and  I0>e  is  the  thermal 
electron  current  to  a  probe  at  plasma  potential  [9].  It  is  assumed  that  the  undisturbed  plasma  is  both  quasineutral  and 
singly-ionized  such  that  the  electron  and  ion  densities  are  equal.  Manipulation  of  Eqs.  (2),  (3),  and  (4)  results  in  the 
fundamental  double  probe  characteristic  given  by  Eq.  (5). 
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The  derivation  of  Eq.  (5)  makes  no  assumptions  about  the  features  of  the  ion  current  collection  mechanism,  and 
thus  is  applicable  for  all  rP/XD  provided  that  the  electron  distribution  function  is  Maxwellian,  such  that  Eq.  (4)  is 
valid.  If  the  ion  saturation  current  is  independent  of  the  applied  bias  potential  (for  /7+=/2+),  then  Eq.  (5)  reduces  to 
the  original  symmetric  double  probe  formulation  of  Johnson  and  Malter  [12]. 


PLASMA 

Il,e  l2,e 
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Ip - ► 

FIG  1.  Electrical  circuit  of  double  Langmuir  probe  in  plasma,  including  voltage  source  meter  with  ion  and 
electron  currents  at  each  electrode. 

It  has  been  shown  that  the  ion  current  collected  by  a  spherical  or  cylindrical  probe  can  be  represented  by  Eq. 
(6),  where  I0  is  the  ion  current  at  the  sheath  edge  given  by  Eq.  (7)  [9,1 1].  The  fit  parameters  a  and  b  are  functions  of 
rP/XD  and  are  given  by  the  expressions  in  Table  II,  which  have  been  reported  to  produce  correlation  coefficients 
greater  than  0.997  over  the  range  3<rPAD<5 0  [11].  For  the  cylindrical  probe  at  values  of  rP/lD< 3,  the  true  b 
parameter  deviates  from  the  value  given  by  Table  I  and  can  be  approximated  as  having  a  constant  value  of  0.5  over 
the  range  0<  rP/kD< 3  [10]. 
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Table  II.  The  fit  parameters  a  and  b  for  3<rP/^D<50  (from  Ref.  1 1). 


Probe  geometry 

a 

B 

cylindrical 

I.I8-0.00080(;>/;.d)135 

0.0684+(0. 722+0. 928 

x  ™ 

spherical 

1.58+(-0.056+0.816  x  />/AD)'0'744 

-0.933+(0.0148+0.119 

x  rP/xDy0-'25 

The  expression  above  for  the  collected  ion  current  can  be  inserted  into  Eq.  (5)  to  yield  the  final  double  probe 
current  characteristic,  which  is  given  by  Eqn.  (8).  To  make  use  of  this  expression  in  deducing  plasma  parameters 
from  experimental  data,  one  must  first  relate  the  potential  of  electrode  1,  Vh  to  a  directly  measurable  quantity  such 
as  Vp.  This  is  accomplished  by  noting  that  the  probe  as  a  whole  floats  such  that  no  net  current  is  drawn  from  the 
plasma,  as  shown  by  Eq.  (9).  This  equation  can  be  solved  implicitly  to  yield  Vj  as  a  function  of  VP  for  a  given  ion 
species,  electron  temperature,  and  predetermined  values  of  the  fit  parameters  a  and  b. 
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III.  Experimental  Apparatus 

A.  Vacuum  Facilities  and  Low  Power  Hall  Effect  Thruster 

The  investigation  was  conducted  at  the  Air  Force  Research  Laboratory  (AFRL)  in  the  Space  Environmental 
Facility  (SPEF),  a  stainless  steel,  spherical  vacuum  chamber  9.1-m  in  diameter.  SPEF  utilizes  six  48”  diffusion 
pumps  with  liquid  nitrogen  cooled  cold  traps.  The  facility  has  demonstrated  xenon  pumping  speeds  greater  than 
300,000  1/s  and  base  pressure  of  lxl 0"6  torr-xenon.  Two  MKS  Instruments  cold  cathode  gauges  (CCG)  were  located 
on  the  walls  as  shown  in  Fig.  2,  and  background  pressure  was  varied  by  changing  the  number  of  operating  diffusion 
pumps. 

A  Velmex  motion  control  system  consisting  of  three-axis  translation  and  a  rotation  stage,  illustrated  in  Fig.  3, 
was  used  to  position  the  plasma  diagnostic  array  through-out  the  Hall  thruster  plume.  The  translation  path  and 
measurement  positions  in  Fig.  3  show  the  measurement  radius  from  0.24  m  to  3.7  m  in  2  degree  increments. 
Following  the  coordinate  system  shown  in  Fig.  3,  this  system  enables  probe  measurements  over  a  5.79-m  X-axis, 
0.5 1-m  Y-axis,  and  4.57-m  Z-axis  range  with  180°  of  rotation.  The  probes  were  mounted  in  parallel  on  a  common 
diagnostic  array,  spaced  approximately  5  inches  apart.  At  each  measurement  location,  the  diagnostic  array  was 
rotated  for  each  probe  such  that  the  probe  was  pointed  toward  the  thruster  center.  Although  this  leads  to  minor 
differences  in  position  between  the  three  probes  at  a  measurement  location,  the  data  is  interpolated  in  Igor  Pro  using 
Delaunay  triangulation.  Probe  positions  were  accurate  to  within  2.5  cm  throughout  each  plume  scan. 

The  ion  source  for  these  experiments  was  a  low-power  laboratory  Hall  thruster.  Prior  to  Langmuir  probe  data 
collection,  the  thruster  was  fired  for  over  one  hour  after  initial  start-up.  Thruster  telemetry  was  monitored  during 
probe  sweeps,  and  exhibited  negligible  deviation  from  steady-state  operation. 
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FIG  2.  Illustration  of  SPEF  with  experimental  setup  within  (left)  and  top-down  view  of  diffusion  pump  and 
CCG  locations  (right). 
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FIG  3.  Illustration  of  SPEF  motion  control  stages  and  axes  (left).  Probe  array  translation  path  and 
measurement  locations  (right). 
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B.  Double  Langmuir  Probes 

The  three  cylindrical,  double  Langmuir  probes  used  in  this  investigation  were  designed  for  OML  analysis,  and 
scaled  to  enable  far- field  measurement  over  a  wide  range  of  number  density  spanning  four  orders  of  magnitude  from 
1016  m"3  in  the  central  plume  to  1013  m"3  in  the  outer  periphery.  Probe  dimensions  are  listed  in  Table  1  and  the  probe 
array  is  shown  in  Fig.  4.  Using  a  single  analysis  method  such  as  OML  required  three  or  more  double  probe 
geometries.  Likewise,  using  any  single  probe  with  multiple  analysis  methods,  including  OML  (rP/A,D<3)  and  thin- 
sheath  (rP/^D>10),  would  limit  the  regions  where  accurate  measurements  can  be  acquired  due  to  low  signal,  poor 
spatial  resolution,  or  blending  of  analysis  methods  where  ?><rP/)iD<\0. 

Cylindrical  probe  electrodes  were  tungsten  rod  and  extended  from  the  diagnostics  array  housing  through  an 
alumina  sleeve.  The  rotation  stage  oriented  each  probe  toward  the  thrust  axis  at  the  Hall  thruster  exit  plane  where 
(X,  Z)  =  (0,  0),  such  that  plasma  flow  was  aligned  with  the  probe  electrode  axis.  Past  studies  have  shown  the  double 
probe  plasma  I-V  characteristic  may  be  analyzed  as  a  stationary  plasma  in  this  orientation,  provided  the  electrode 
length  is  much  greater  than  the  diameter  [13].  In  this  study,  all  electrode  length  to  diameter  ratios  exceeded  10:1, 
and  are  sufficient  to  meet  this  criteria.  A  Keithley  2410  sourcemeter  was  swept  from  -15  V  to  +15  V  using  the 
circuit  schematic  illustrated  in  Fig.  1. 


Table  1  Double  Langmuir  probe  dimensions 


Double  Langmuir  Probe 

Electrode  Diameter  [mm] 

Electrode  Length  [cm] 

Electrode  Spacing  [cm] 

Probe  1 

0.254 

5.08 

4.45 

Probe  2 

3.175 

12.2 

5.72 

Probe  3 

9.525 

19.70 

5.72 

FIG  4.  Photographs  of  the  diagnostics  array  (left)  with  close-up  double  Langmuir  probe  2  (right). 
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IV.  Results  and  Discussion 

Double  Langmuir  probe  results  were  normalized  to  a  consistent  value  and  plotted  with  respect  to  thruster 
diameter,  ranging  150  diameters  downstream  and  to  80  diameters  on  the  periphery.  Probe  1  was  designed  for  the 
highest  density  regions,  and  exhibited  minimal  signal  to  noise  beyond  approximately  40  thruster  diameters.  Probe  2 
was  designed  to  span  the  entire  range  of  plasma  properties  throughout  the  plume  with  OML  analysis.  Probe  3  was 
designed  for  the  low  density  regions  on  the  plume  periphery,  and  due  to  its  size  suffered  from  poor  spatial  resolution 
within  20  thruster  diameters. 

Number  density  calculated  using  the  analysis  technique  described  in  Ref.  7  is  compared  to  OML  in  Fig.  5  for 
Probes  2  and  3.  This  qualitative  assessment  indicates  the  number  densities  calculated  with  Ref.  7  are  consistent 
between  Fig.  5a  and  Fig.  5c,  corresponding  to  Probe  3  and  2  respectively.  For  both  probes,  number  density  in  the 
central  plume  in  Figs  5a  and  5c  calculated  with  Ref.  7  extend  further  downstream  compared  to  OML  analysis  in 
Figs.  5b  and  5d.  Comparison  the  central  plume  in  Figs.  5c  and  5d  indicates  Probe  2  has  better  agreement  between 
Ref.  7  and  OML,  whereas  Figs.  5a  and  5b  show  Probe  3  results  in  similar  features  between  Ref.  7  and  OML  on  the 
periphery.  These  trends  in  OML  analysis  are  consistent  with  expected  results  due  to  the  poor  spatial  resolution  of 
Probe  3  nearest  the  thruster  and  the  reduced  signal  of  Probe  2  on  the  far  wings  of  the  plume.  The  agreement 
between  Probe  2  and  3  using  the  analysis  of  Ref.  7  indicates  it  accounts  for  sheath  physics  across  a  range  of  rP/XD , 
and  is  suitable  for  the  transitional  regime  where  OML  is  not  valid,  when  rPAD> 3. 


X  [Thmster  Diameters] 

(d) 


Normalized  Ion 
Number  Density 


Normalized  Ion 
Number  Density 


FIG  5.  Comparison  of  normalized  number  density  of  double  Langmuir  Probe  3  using  analysis  of  (a)  Ref.  7 
and  (b)  OML,  and  of  Probe  2  using  analysis  of  (c)  Ref.  7and  (d)  OML  at  low  background  pressure. 
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In  Fig.  6,  number  density  calculated  using  the  analysis  technique  described  in  Ref.  7  is  compared  to  OML  for 
Probes  1  and  2,  where  the  signal  to  noise  of  Probe  1  is  adequate  for  analysis  of  plasma  properties.  This  qualitative 
assessment  indicates  the  number  densities  calculated  with  Ref.  7  are  consistent  between  Fig.  6a  and  Fig.  6c, 
corresponding  to  Probe  2  and  1  respectively.  Comparison  the  central  plume  in  Figs.  6c  and  6d  indicates  Probe  1  has 
better  agreement  between  Ref.  7  and  OML,  whereas  Figs.  6a  and  6b  show  Probe  2  results  in  conflicting  features 
between  Ref.  7  and  OML  in  the  highest  density  region  nearest  the  thruster.  These  trends  indicate  the  Probe  2  OML 
analysis  in  Fig.  6b  are  influenced  by  reduced  spatial  resolution  of  Probe  2  nearest  the  thruster  and  that  the  probe  is 
no  longer  in  the  OML  analysis  regime.  The  agreement  between  Probe  1  and  2  using  the  analysis  of  Ref.  7  indicates 
it  accounts  for  sheath  physics  across  a  range  of  rP/lD,  and  is  suitable  for  the  transitional  regime  where  OML  is  not 
valid,  when  rP/XD> 3. 

Centerline  number  density  from  the  three  double  probe  configurations  are  calculated  with  thin  sheath, 
parametric  analysis  in  Ref.  7,  and  OML  in  Figs.  7  to  9.  Trends  indicate  the  analysis  of  Ref.  7  is  in  agreement  with 
thin  sheath  analysis  out  to  approximately  25  to  45  thruster  diameters,  and  in  agreement  with  OML  analysis  beyond 
approximately  50  diameters.  This  suggests  the  analysis  of  Ref.  7  is  bridging  the  plasma  properties  from  OML 
analysis  to  thin  sheath,  and  is  advantageous  for  determination  of  plasma  properties  across  wide  regions  of  the  Hall 
thruster  plume.  Based  on  past  simulations  of  the  plume,  the  location  between  approximately  30  to  50  thruster 
diameters  downstream  of  the  exit  plane  corresponds  to  the  transition  region  where  3<fp/ad<\  0. 


X  [Thmster  Diameters] 


X  [Thmster  Diameters] 


(c) 


(d) 


FIG  6.  Comparison  of  normalized  ion  number  density  of  double  Langmuir  Probe  2  using  analysis  of  (a)  Ref. 
7  and  (b)  OML,  and  of  Probe  1  using  analysis  of  (c)  Ref.  7  and  (d)  OML  at  low  background  pressure. 
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Thruster  Diameters  [  -  ] 


FIG  7.  Comparison  of  normalized  ion  number  density  of  double  Langmuir  Probe  1  using  thin  sheath,  Ref.  7, 
and  OML  at  low  background  pressure. 


Thruster  Diameters  [  -  ] 

FIG  8.  Comparison  of  normalized  ion  number  density  of  double  Langmuir  Probe  2  using  thin  sheath,  Ref.  7, 
and  OML  at  low  background  pressure. 
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Thruster  Diameters  [  -  ] 


FIG  9.  Comparison  of  normalized  ion  number  density  of  double  Langmuir  Probe  3  using  thin  sheath,  Ref.  7, 
and  OML  at  low  background  pressure. 
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The  centerline  number  densities  all  three  probes  calculated  using  Ref.  7  in  Figs.  7  to  9  are  shown  in  Fig.  10. 
This  shows  good  centerline  agreement  between  all  three  probes  to  40  thruster  diameters  downstream,  and  highlights 
issues  in  Probe  1  due  to  low  signal  to  noise.  Consistent  values  between  Probe  2  and  Probe  3  are  further  indication 
that  the  analysis  in  Ref.  7  accounts  for  probe  sheath  effects.  Fig.  11  shows  the  electron  temperature  on  thruster 
centerline  for  all  three  probes  corresponding  to  data  in  Fig.  10.  While  there  is  considerable  uncertainty  in  electron 
temperature  calculations,  it  shows  consistent  agreement  between  Probes  2  and  3  from  20  to  150  thruster  diameters, 
is  consistent  with  past  studies  using  this  low-power  Hall  thruster. 

The  agreement  between  Probes  2  and  3  on  thruster  centerline  in  Fig.  10  are  further  evaluated  in  Figs.  12  to  17 
for  low  and  high  facility  background  pressures  on  centerline,  0=30°,  and  0=60°  from  the  thruster  centerline  axis. 
Comparisons  between  Ref.  7  and  OML  reveal  significantly  better  agreement  between  parameters  calculated  with 
Ref.  7.  The  exception  is  for  Probe  3  within  40  thruster  diameters  when  0=30°,  which  is  attributed  to  poor  spatial 
resolution.  Beyond  40  thruster  diameters,  calculations  using  Ref.  7  and  OML  are  in  good  agreement  for  0=30°  and 
0=60°  for  both  low  and  high  pressure  cases. 

Comparisons  of  calculated  plasma  properties  using  OML,  thin  sheath,  and  Ref.  7  do  not  necessarily  reduce 
measurement  uncertainty,  but  provide  confidence  in  measurements  with  a  fixed  probe  design  using  a  single  analysis 
technique  in  Ref.  7.  The  wide  range  of  plasma  properties  over  this  expansive  range  of  the  far- field  Hall  thruster 
plume  is  highly  advantageous  for  investigations  of  a  Hall  thruster  or  ion  thruster  plume,  due  to  the  increased  range 
of  plasma  properties  that  can  be  measured,  the  reduced  complexity  of  analysis,  and  the  flexibility  to  study  a  wide 
range  of  properties  in  a  time-varying  plasma  with  variation  in  sheath  characteristics. 


Thmster  Diameters  [  -  ] 

FIG  10.  Comparison  of  normalized  ion  number  density  of  double  Langmuir  Probes  1,  2,  and  3  at  0=0°  using 
analysis  of  Ref.  7  at  low  background  pressure. 


Thruster  Diameters  [  -  ] 

FIG  11.  Normalized  electron  temperature  of  double  Langmuir  Probes  1,  2,  and  3  at  0=0°  using  analysis  of 
Ref.  7  at  low  background  pressure. 
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Thruster  Diameters  [  -  ] 

FIG  12.  Normalized  ion  number  density  of  double  Langmuir  Probes  2  and  3  at  0=0°  using  analysis  of  Ref.  7 
and  OML  at  low  background  pressure. 


0  20  40  60  80  100  120  140  160 

Thruster  Diameters  [  -  ] 

FIG  13.  Normalized  ion  number  density  of  double  Langmuir  Probes  2  and  3  at  0=30°  using  analysis  of  Ref.  7 
and  OML  at  low  background  pressure. 


0  20  40  60  80  100  120  140  160 

Thruster  Diameters  [  -  ] 


FIG  14.  Normalized  ion  number  density  of  double  Langmuir  Probes  2  and  3  at  0=60°  using  analysis  of  Ref.  7 
and  OML  at  low  background  pressure. 
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Thruster  Diameters  [  -  ] 

FIG  15.  Normalized  ion  number  density  of  double  Langmuir  Probes  2  and  3  at  0=0°  using  analysis  of  Ref.  7 
and  OML  at  high  background  pressure. 


Thruster  Diameters  [  -  ] 


FIG  16.  Normalized  ion  number  density  of  double  Langmuir  Probes  2  and  3  at  0=30°  using  analysis  of  Ref.  7 
and  OML  at  high  background  pressure. 
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FIG  17.  Normalized  ion  number  density  of  double  Langmuir  Probes  2  and  3  at  0=60°  using  analysis  of  Ref.  7 
and  OML  at  high  background  pressure. 
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V.  Summary  and  Conclusions 

A  new  double  Langmuir  probe  analysis  technique  based  on  fits  to  Laframboise’s  data,  and  self-consistently 
accounts  for  probe  sheath  expansion  in  the  local  plasma.  The  technique  extends  past  parametric  fits  of  single 
Langmuir  probe  data  to  the  double  probe  configuration  in  a  low-temperature  plasma.  The  methods  enable  double 
probe  measurement  of  electron  temperature  and  plasma  density  over  a  wide  range  of  densities  without  advance 
knowledge  of  the  probe  radius  relative  to  Debye  length.  This  has  implications  for  laboratory  experiments  and  on- 
orbit  measurements. 

The  technique  was  evaluated  with  plume  measurements  of  a  low-power  Hall  thruster  using  three  double 
Langmuir  probe  designs.  Plasma  properties  calculated  with  the  new  analysis  technique  were  consistent  with  OML 
and  thin  sheath  analysis,  and  demonstrated  agreement  between  the  three  probes  throughout  the  plume.  This  was 
successfully  verified  at  two  facility  background  pressures  and  number  densities  spanning  four  orders  of  magnitude 
out  to  150  thruster  diameters  from  the  thruster  exit  plane.  The  study  suggests  this  double  Langmuir  probe  technique 
is  well-suited  for  low-temperature  plasmas  in  the  far-field  electric  propulsion  plume,  and  has  several  key  advantages 
over  traditional  analysis  techniques,  including:  the  ability  to  measure  a  wide  range  of  plasma  density  with  a  fixed 
probe  design  and  a  single  analysis  technique,  reduced  complexity  in  data  analysis  over  a  wide  range  of  ry^,  and  the 
flexibility  to  study  a  wide  range  of  properties  in  a  time-varying  plasma  with  variation  in  sheath  characteristics. 
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Background  and  Motivation 


Background 

•  Double  probe  theory  based  on  theoretical  analysis  at  AFRL1 

•  Valid  for  low  temperature  plasma 

•  Based  on  analytical  fits  to  Laframboise’s  results,  extended  to  double  probe  geometry 

Advantages  over  conventional  techniques 

•  Self-Consistent  — >  Accounts  for  sheath  expansion 

•  Simple  — ►  Enables  single  analysis  method  with  fixed  probe  design 

•  Enabling  — ►  Expands  measurement  range  for  transient  plasmas 


[1]  Beal,  B.  E.,  Johnson,  L.,  Brown,  D.  L.,  Blakely,  J.  M.,  Bromaghim,  D.,  “Improved  Analysis  Techniques  for  Cylindrical 
and  Spherical  Double  Probes,”  (Accepted  by  Rev.  Sci.  Instrum.  2012). 
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w  Probe  Theory  (1/3) 

Standard  Langmuir  Probe  Analysis  Techniques 


Orbital  Motion  Limited  (OML) 

•  Valid  for  small  diameter  probes,  low  density  plasma 

•  Plasma  sheath  infinitely  large  compared  to  the 
probe  dimensions  (Tp/AD<3) 

•  Sheath  does  not  limit  the  penetration  of  the  electric 
field  into  the  bulk  plasma 

•  Probe  collects  all  ions  whose  momentum  relative  to 
the  probe  surface  is  insufficient  to  escape  the 
electric  field 

Thin  Sheath  Analysis 

•  Valid  for  large  diameter  probes,  high  density  plasma 

•  Plasma  sheath  small  relative  to  the  probe 
dimensions  (Tp/AD>10) 

•  Probe  collects  all  ions  that  enter  the  sheath 


*  J.G.  Laframboise,  Theory  of  Spherical  and  Cylindrical 
Langmuir  Probes  in  a  Collisionless,  Maxwellian 
Plasma  at  Rest,  (University  of  Toronto  Institute  for 
Aerospace  Studies,  Toronto,  1966),  Report  No.  100. 
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Interpretation  requires  knowledge  of  the  relation  between  the  ion  current 
collected  by  a  biased  electrode  and  the  local  plasma  parameters 
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w  Probe  Theory  (2/3) 

Standard  Langmuir  Probe  Analysis  Techniques 


Laframboise’s  Results  for  Collisionless 
Plasma2 

•  Several  analytical  fits  developed,  valid  over  a  limited  range  of 

rp/AD 

•  Iterative  approach  to  determine  properties 

•  Requires  prior  knowledge/assumptions  about  local  plasma 

•  Multiple  regimes  and  transition  complicates  analysis 

•  Typically  achieve  -50%  uncertainty 

•  Analytical  parameterization  for  single  electrode  probes  that  is 
applicable  over  a  wide  range  of  rp/AD  for  cold  ions  (T/Te«1)3’4’5 

•  Encompasses  wide  range  of  plasma  properties  (OML,  thin 
sheath,  and  transition) 

Extend  Parameterization  to  Double 

Probe  Geometry 


[2]  J.G.  Laframboise,  Theory  of  Spherical  and  Cylindrical  Langmuir  Probes  in  a  Collisionless,  Maxwellian 
Plasma  at  Rest,  (University  of  Toronto  Institute  for  Aerospace  Studies,  Toronto,  1966),  Report  No.  100. 

[3]  A.  Karamcheti  and  Ch.  Steinbruchel,  J.  Vac.  Sci.  Technol.  A,  17,  3051  (1999). 

[4]  Ch.  Steinbruchel,  J.  Vac.  Sci.  Technol.  A,  8,  1663  (1990). 

[5]  G.  Narasimhan  and  Ch.  Steinbruchel,  J.  Vac.  Sci.  Technol.  A,  19,  376  (2001). 


[Ref  1] 


FIGURE  21  ION  CURRENT  4.  VS  PROBE  POTENTIAL  FOR  VARIOUS  RATIOS 
OF  PROBE  RADIUS  TO  ELECTRON  DEBYE  LENGTH;  ION- 
ATTRACTING  SPHERICAL  PROBE;  T*/T.  •  0;  ELECTRONS 
REFLECTED  BY  PROBE  SURFACE;  OBTAINED  FROM  NUM¬ 
ERICAL  SOLUTION  OK  THE  ALLEN.  BOYD.  AND  REYNOLDS 
EQUATION.  DOTTED  CURVE  SHOWS  TRAPPED-ORBIT 
BOUNDARY. 
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w  Probe  Theory  (3/3) 

Improved  Double  Probe  Technique 


Extend  Parameterization  Technique  to  Double  Probe 
Geometry1 

•  Self-consistent  manner  that  does  not  require  knowledge  of  the  local  plasma  potential 
(unlike  the  single  probe) 

•  Find  electron  temperature  and  plasma  density  over  wide  range  of  rp/AD  for  cold  ions 

•  Double  probe  circuit  with  Kirchoff’s  Laws 

•  Accounts  for  finite,  but  non-negligible  extend  of  plasma  sheath  in  transition  region 


PLASMA 


(Kiithk-y  2410) 
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electron  distribution  function  is  Maxwellian 


Fitting  Parameters  for  Analytical  Parameterization3’4’5 
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Experimental  Overview  (1/2) 

Facilities  and  Plasma  Environment 


SPEF  vacuum  chamber 

—  9.1-m  diameter  stainless  steel 
sphere 

—  Six  48”  diffusion  pumps,  greater  than 
300,000  L/s  on  Xe 

—  Velmex  Motion  control  system  w/  3- 
axis  translation,  180  degree  rotation 


axis 


57 


m 


•  Low  power  Hall  thruster  ion  source 

•  Two  background  pressures  (vary  #  DPs) 

•  Plasma  conditions  expected  to  range  from 
OML  to  thin  sheath  based  on  past 
measurements  with  this  thruster 
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Langmuir  Probe  Design 

—  3  Probes  scaled  to  measure  for 
measurement  throughout  plume 

—  Evaluate  w/  new  technique  to  compare  with 
thin  sheath  and  OML  (where  valid) 

—  Tungsten  rod  electrodes,  alumina  tube 

—  Keithley  2410  sourcemeter  used  to  sweep 
electrodes  from  -15  V  to  +15  V 

—  Measure  150  thruster  diameters 
downstream  (+/-  80  diameters  on  periphery) 
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Double  Langmuir 
Probe 


Probe  1 


Probe  2 


Probe  3 
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Langmuir  Probe  Plume  Data 

Plume  Comparison  of  Probes  2  and  3, 1.3X1  O'6  torr-Xe 


Qualitative  Comparison 

•  Probe  2  and  probe  3  in  agreement 
using  new  analysis 

•  Probe  3  shows  better  agreement 
on  periphery 

•  Probe  2  shows  better  agreement 
in  central  plume 

•  Poor  agreement  using  OML  theory 
reveals  large  uncertainty  typical  of 
Langmuir  probe  measurements 
(attributed  to  sheath  physics  with 
local  plasma) 
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Langmuir  Probe  Plume  Data 

Plume  Comparison  of  Probes  1  and  2, 1.3X1  O'6  torr-Xe 


Qualitative  Comparison 

•  Probe  1  and  probe  2  in  agreement 
using  new  analysis 

•  Probe  2  shows  poor  agreement 
throughout  near-field 

•  Probe  1  shows  better  agreement 
in  central  plume  using  OML 

•  Poor  agreement  using  OML  theory 
reveals  large  uncertainty  typical  of 
Langmuir  probe  measurements 
(attributed  to  sheath  physics  with 
local  plasma) 
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Langmuir  Probe  Plume  Data 

Comparison  of  Analysis  Techniques,  1.3X1  O'6  torr-Xe 


New  analysis  consistent  with 
thin  sheath  in  near-field  and 
OML  in  far-field 

Probe  1  not  valid  in  far-field  beyond  40 
diameters  (poor  signal  to  noise  for  small 
diameter  probe) 

Probe  3  not  valid  in  near  field  due  to  poor 
spatial  resolution  (large  probe  with  length 
several  times  thruster  diameter) 

Data  shows  that  a  single  probe  is  not  adequate 
to  accurately  measure  full  range  of  properties 
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w  Langmuir  Probe  Plume  Data 

Assessment  of  Technique  wrt  Probe  Geometry,  1 .3X10'6  torr-Xe 
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Data  shows  plume  scattering 
and  ion  migration  at 
approximately  30  degrees  off- 
axis  is  accounted  for  with  new 
method  (attributed  to 
accounting  for  sheath  effects) 

Demonstrates  consistent 
agreement  between  Probe  2 
and  3  throughout  the  plume 


NOTE:  Probe  1  did  not  have 
adequate  signal  to  noise  for 
comparison  at  30,  60  degrees 
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High  pressure  shows  new 
analysis  accounts  for  sheath 
effects,  consistent  with  low 
pressure  data 

New  analysis  consistently 
shows  agreement  in  number 
density  within  10% 

Future  efforts  may  focus  on 
uncertainty  analysis  to  reduce 
50%  error  bars 
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Summary  and  Conclusions 


New  double  probe  analysis  technique  evaluated  with 
measurement  of  Hall  thruster  plume  using  3  probe  designs 


Comparison  indicates  new  analysis  technique  successfully  accounts  for 
probe  sheath  physics  over  wide  range  of  rpIK0  throughout  the  plume 

Probes  show  better  agreement  throughout  the  plume  with  new  analysis 

Analysis  advantageous  for  transition  regions  between  thin  sheath  and  OML 
analyses 

Investigation  supports  using  a  fixed  probe  design  with  single 
analysis  over  a  wide  range  of  rp/AD  in  low  temperature  plasma 
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